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The problem of random motion of ID quantum reactive harmonic oscillator (QRHO) is formulated in 
terms of a wave functional regarded as a complex probability process ^stc (x, t; {£}) in an extended space 
S = R 1 (8 R{£}- In the complex stochastic differential equation (SDE) for ty s tc (x, t; {£}) the variables 
are separated with the help of the Langevin-type model SDE introduced in the functional space R{^y- The 
complete positive Fokker-Plank measure of the space R{^y is obtained. The average wave function of roaming 
QRHO is obtained by means of functional integration over the process ^stc {x, t; {£}) with the complete 
(f) ' Fokker-Plank measure in the space R{£}- The local and averaged transition matrices of roaming QRHO 

were constructed. The thermodynamics of nonrelativistic vacuum is investigated in detail and expressions 
OS ■ for the internal energy, Helmholtzian energy and entropy are obtained. The oscillator's ground state energy, 

its shift and broadening are calculated. 
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1 Introduction 



From the point of view of classical dynamics all phenomena, that are described in the framework of conventional 
theory of quantum mechanics are stochastic processes. The natural symmetry between the Schrodinger and 
Fokker-Plank equations was used both for the formulation of quantum mechanics as a classical stochastic theory 
pi] and quantization of classical theory by introduction of the concept of random process [0 . 

Note, that in all representations the main object of quantum mechanics, the wave function, is deterministic. 
It is worthwhile to stress that the deterministic nature of a physical theory is a consequence of the symmetry 
of basic equations with respect to time inversion || . 
, It is known that the observability is an essential property of objects that is realized when the contact with 

an external world is established. Therefore, all objects under study are open physical systems. This conclusion 
implies, in particular, that it is impossible to conform in any reasonable way the probabilistic, irreversible nature 
of measurements to the reversibility of quantum theory concepts H . 

Moreover, the stochastic nature is the basic property of the physical world. This assertion becomes especially 
conclusive if one allows for the real existence of physical vacuum in the nature that is the cause of casual 
quantum jumps in the isolated systems ||. In other words, all objects that are investigated taking into account 
the interaction with fluctuations of the physical vacuum are open physical systems. Whereas the isolated 
systems are studied by means of methods of conventional quantum mechanics, the behaviour of open systems, 
in particular of systems that are in equilibrium with the environment, is usually described by the laws of 
thermodynamics. 

Since due to permanent fluctuations the open system has no any definite quantum state, i.e., the state vector 
of the system randomly changes with time, the construction of a theory alternative to the reversible quantum 
mechanics becomes urgent. One can construct such a theory if, in particular, the random processes taking 
place in the system are directly included in the foundations of the theory and are not regarded as external 
perturbations. 

In the present work the problem of stochastic quantum mechanics is investigated using the model of random 
motion of ID quantum reactive oscillator (QRHO) . Mathematically the problem of an evolution of quantum 
states is formulated in the framework of probabilistic complex process in an extended space with stochastic 
space-time continuum. For roaming classical oscillator the conditions of reducibility of a model stochastic 
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differential equation (SDE) to nonlinear Langevin SDE are investigated. With the help of model SDE the 
variables in the SDE for a wave functional - the complex probability process, have been separated. Based on 
the Langevin-type SDE the Fokker-Plank equation is derived and the measure of functional space is obtained. 
For roaming QRHO the average wave function is constructed in the form of functional integral with positive 
measure. The local stochastic and averaged transition matrix of the roaming QRHO is calculated analytically. 
The probability of " vacuum- vacuum" transition is calculated numerically and its behavior is shown versus the 
fluctuation parameter. The thermodynamics of nonrclativistic vacuum is builded and the expressions for the 
main thermodynamic potentials are obtained. 



2 Formulation of the problem 

Let us consider a closed system "quantum object -(-thermostat" and assume that the quantum object moves in 
the Euclidean space R 1 and the state of the thermostat is characterized by an unlimited set of modes in the 
functional space R{£}, where {£} = £ (i, {W}) is a complex functional of some random process {W} = W (t). 
One can mathematically describe the wave state of the system with the help of a functional - a complex 
probabilistic process ^ s tc (x, t; {£}), given on the extended space S = R 1 (g> R{^}- The equation of evolution for 
the wave functional can be written in the most general case as: 

idt^stc = H (x,t;{W})^ stc , -co < x, t < +co, (1) 



where dt = d/ dt is the total derivative of the random functional (see (|10[)) with respect to time, and 
H (x, t; {W}) is respectively the stochastic operator of the system evolution. 

Below, as a specific example, we shall consider the case of roaming ID QRHO, the operator of evolution of 
which has the following form: 

H(x,t;{W}) = ~[-d 2 x + n 2 (t;{W})x 2 ], d 2 x ^d 2 /dx 2 . (2) 
Let the frequency be given as: 

{l(t;{W}) = n (t) + Sl 1 (t;{W}), (3) 
where the regular part Qq (t) > satisfies the following boundary conditions: 

lim n (t) = n in(out) > 0. (4) 

Note, that when f^i (t; IW"}) = 0, the equations (Q)-^) describes the case of regularly moving QRHO with 
asymptotic wave function (0), that admits exact solution (see 0). In the expression (||) fli (t; {VF}) respectively 
denotes the random part of frequency and satisfies the following conditions: 

lim fii (t; {W}) = 0, lim Q x (t; {W}) ^ 0. (5) 

t — > — 00 t — >+oo 

Next, assume that the total wave functional meets the natural boundary conditions: 

lim * stc (x, t; U» = lim d x V sfc (x, t; {£}) = 0. (6) 

— >oo \x\ — >oo 

We denote as (n\x, t; {£}) the total wave functional that is evolved in the Euclidean subspace R\ n from 
the purely n-th vibrational state of the quantum harmonic oscillator: 

-,(+) 



^JJ(n\x,t;{^}) - V in (n\x,t) 
t — *— 00 



1 -| 1/2 

2™n! 



exp |-i + n in t - ifi m a; 2 | H n (^/D~x^j , (7) 



0,1,2. 
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The objectives of the present investigation are: 

a) to establish the criteria for separation of variables in SDE and determine the random wave 
functional (n\x, t; {£}) in an explicit form; 

b) to calculate the wave function of roaming QRHO 

*£ ) (n\x,t) = (¥£(n\x,t ] {0)) m , (8) 

where the brackets (. . .) r^i stand for the functional integration including the integration over the distribution 
of coordinate £ in the functional space R{£} at instant t; 

c) to calculate the transition matrix S^" m of the roaming QRHO. 

3 Derivation of SDE for roaming classical oscillator 

We begin our consideration with a second-order equation 

e + n 2 (i;{W/})£ = 0, (9) 

where 

£ = dtt (t; {W}) , d£ (t; {W}) = (dtt + ±5^A dt + (6 w t) dW (t) . (10) 

Remember that the second formula in ([l0|) denotes the Ro differential of the random functional £ (t; {W 7 }) 
(see 0), where the functional derivative is determined in the usual way: 

Taking into account the conditions we find the asymptotical solution of in the (in) channel to be 

Ht;{W}) ~ exp(tfW). (12) 

t — > — OQ 

One is to stress here, that the equation (^|) makes sense and describes the behaviour of roaming classical 
oscillator only when it is reducible to the Langevin-type SDE. 

Theorem: There exists a set of nonsingular junctionals £ (t; {W}) that reduce the equation to a nonlinear 
complex Langevin-type SDE, or, that is the same, to a system of two real nonlinear SDE. 

Proof: Assume that the solution of model equation (Q) has the following form: 

£ (t; {W}) = £ (*) exp ( J <p (*'; {W 1 }) dA , (13) 

where £o (t) is the solution of (|^) with regular frequency f2o(t) (see Q). Consider the first total derivative 
of the random functional (|l^) with respect to time. Taking into account ( |l0| ) one has 

dti (*; {W}) = [Cot (t)/ £o (t) + (t; {^}) + 
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- DC 



0(* , ;W)A , + (dt^(*)) / <W(*';{W}H*' 

— OO 

£ot (t) = d t £o (t) . 



t(t;{w}), (w) 



Since for all t' < t 

W(t';{W})| t , <t = 0, (15) 
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we can rewrite the expression ( Jl4| ) as 

d t (; (t; {W}) = [Cot (*)/ £o (t) + (*; {W})] £ (t; {W}) = dt£ (t; {W}) 
Substituting the solution (|l^) to (JsJ> and taking into account ([l6]) one can find: 

+ 2 [fo (*)/ Co (*)] 4> + 4> 2 + nl (t) + U (t; {W}) = 0, 
where the notation was made 

U (t; {W}) = fi? (i; {W}) + 2fii (t; {W}) fi (f) . 
In particular, from (|l5|) and equation (|lj]) the following relations results 

8wt(t;{W}) = 0, 8w{dtt(t;{W})}?0 

from which it is seen that the operators Sw and dt do not commute. 
Making the substitution 

ct>{t;{W})=^{t;{W})-^t (t)/&(*). 

one obtain finally from jl7| ) 

$ + 3> 2 + ft 2 , (*) + + F (t; {W}) = 0. 
In SDE ( |2l| ) the following notations were made 

U (t; {W}) = F + F (i; {W}) , F = (U (t; {W})) , 



(16) 
(17) 
(18) 
(19) 

(20) 
(21) 
(22) 

(23) 
(24) 

After substitution of ( p4f ) to (|2l]) and separation of real and imaginary parts we have the following system 
of SDE: 



and in 122) the averaging of expression (18) is made over the ensemble 



The solution (|21[) of SDE, taking into account (|12[), satisfies the initial condition of the type 

lim $ (t; {W}) = ifl in , 

t — * — oo 

that, in its turn, suggests the complex solution 

$ (t; {W}) = 9 (t; {W}) + tip (t; {W}) . 



6 + 9 2 - if 2 + U (t) + F (i; {W}) = 0, U (t) = n 2 (t) + F > 0, 



(^ + 2(^61 = 0. 



The system of SDE (25)-(E6J) can be rewritten in a vector equation: 



$+K(<;{*})+F(i;{iy}) = 0, 
where the vectors have the following projections 

* (t; {W}) = {6, V } , K (t; {*}) = { [9 2 - <p 2 + U (i) + F ] ; 26} , 
F(i;{^}) = {F(t;{iy});0}. 

The theorem is proved. 

For subsequent investigations of SDE (|l])-(||) it is convenient to write SDE (JsJ) in the form 



(25) 
(26) 

(27) 



(28) 



£ (t; {W}) = a (i; {W}) exp [ir (f; {W})] 



(29) 
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where a (i; {W}) and r (t; {W}) are yet unknown real functionals: 



Im cr (<; {W}) =Imr (i; {W }) = 0. (30) 
Taking into account ([l2]) one can establish the asymptotic behaviour of these functionals: 

t Hrn^ a (i; {W}) = 1, r (t; {W}) Q in t. (31) 

After substitution of ( |29| ) to ([n]) and taking into account (|4|), (|2^) and ([50|) we find that 

6 w a (t; {W}) = 0, 5 w [5f<7(t;{W})] = 0, (32) 

as well as that 

<Vr(t;{W})=0, ^[5 t r(t;{W})] = 0. (33) 
It must be noted that the conditions (|32[)-(p3[) are highly important for our further discussions. 



4 Solution of SDE for complex probabilistic process - the wave func- 
tional 



Now pass to the solution of SDE (n 
Theorem: If the relations (p. 



form of orthonormal in Li (i? 1 ® R{£}) random complex functionals 

i; ,-|l/2 



using the Langevin-type model nonlinear SDE for separation of variables. 
j-(^3) are fulfilled, then the complex SDE has an exact solution in the 



(fW7T) 1/2 



2 n n\a 



vx\> <; / j ,-/ + - | fl, n T- 



'a)' 



(34) 



where r(t;{W}) = r (t; {W})/Cl in , a t = d t a (t; {W}) and r t = d t r(t;{W}). Moreover, L 2 (R 1 ® R { ^ } ) 
denotes the space of square integrable functions on the real axis R 1 with values in some functional space R{£} ■ 
Proof: Changing the variables in 



x y = 



a(t;{W}) 



(35) 



and taking into account (32)-([33|) the equation for complex stochastic process will be as follows: 
id t V stc (y, t; {£}) = U^ySy - + ytt 2 (i; {W}) y 2 | ^ stc (y, t; {£}) , 

Now, let us rewrite the solution of equation ([36j) as 

- , r , % f exp \i2A (t; {W}) y 2 ] 1 1/2 , , r 

Making the substitution of (|37]) into (^6|) and the transformation 

r(t;{W}) 



t -> r 



(36) 



(37) 



(38) 



with regard to (|32|)-(|33]) one finds 
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= ~\ {S 2 y - a 2 [2A - iata^A + 4a" 2 A 2 + ct 2 ^ 2 (f; {W})] y 2 } x , 



(39) 



where the functionals a (t; {W}), r (t; {W}) and A (t; {W}) are still unknown. For their determination we 
shall suppose that the following relations were observed: 



n (t;{W}) = 



° 2 (t;{W})" 



A(t;{W}) = -a t (t;{W})a{t-{W}), 



2A - iata^A + 4a- 2 A 2 + a 2 Q 2 (i; {W}) 



1 



Hi 



(40) 
(41) 
(42) 



It is easy to see that the system of equations ([tt|)-([42|), taking into account (|28|) and relations fl32|)-(|33D, 
is equivalent to the model equation (^) and, thereby to Langevin-type nonlinear SDE ( ^l| ) with the complex 
initial condition (E3h. Using the expressions (M0)-(M3) one can finally obtain from (By): 
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iSrX (y, t) = - (s 2 y + n 2 n y 2 ) x (y, t) . 



(43) 



The equation ( |43| ) describes an autonomous system of the quantum harmonic oscillator on the stochastic 
space-time continuum. The solution of equation (j43|) is obtained by conventional methods and has the form 



(n m /n) 1/2 



2 n nl 



1/2 



exp j-i fn + Q in T - i^ in y 2 | H n (^\/ih^yj 



(44) 



Combining (|37|) and (|44|) for the wave functional, evolving from the initial pure state ((7J) , it is easy to find the 
expression (|3~i|). Taking into account (^)-(^) one can find the relation between the stochastic time r (t; {W}) 
and the apparent time, the natural parameter t, 



dt' 



CT 2 (t';{W'}) 

— OO 

Now take note of the very important property of the wave functional 
y+}(m\x,t; (n\x, t; = S mn , 



+oc 



dx, 



pointing to the fact that the basis formed in the space Li (i? 1 ® R{£}) is orthonormal. 
The theorem is proved. 



(45) 



(46) 



5 Derivation of Fokker-Plank equation for conditional probability 

Let us consider the functional of the form: 

P(<M*',t') = (6 [*(t)-* (*')]>, (47) 

where = * (t; {W}) is the solution of SDE (^?j). After the differentiation of functional (fflj ) over the 

time and use of (^?j) one can obtain: 



G 



d t p (*, t|#', <') = -a* (* 5 [* (*) - * (*')]) = 



9* {K (t- {*}) P (*, t\#J) + (F (t; {W}) 5 [* (t) - * (*')]>} . (48) 



To obtain the equation for conditional probability MTH in an explicit form one has to specify the stochastic 
frequency F(t;{W}). Consider the most common case when the stochastic component ( |2§| ) is the Gaussian 
random function F (t; {W}) = W = F (t). This assumption implies that one can completely define F (t) with 
the help of correlation function that is a white noise type correlator in case when F (t) changes faster than the 
solution 9 (t) 

(F (t) F (£')) = 2eS (i — t') , (F(t)) = Q, e > 0. (49) 
Now, using the Vick's theorem (see j?j]) 

{F(t)N(t;{F}))=^(^) p > (50) 
where N (t; {F}) is an arbitrary functional of F (t), one can write the following expression 

(F (t; {W}) 5 [$ (t) - * (*')]) = 2d e ■ 5 [* (t) - * (*')]) • (51) 

Due to the stochasticity of 9 it) the variational derivative of F (t) is e • sgn(t — £') + 0(t — t'). After 
common regularization procedure (in the sense of Fourier decomposition) one can find the value at t = if: 
e ■ sgn(0) = is. Taking into account the aforesaid we finally obtain the following Fokker-Plank equation for 
conditional probability: 

2 

d t P(*,t\&,t')=J2 d** [ K i (*;{*» + £yd*JP (<&,*!*', i'), (52) 

where 

en = £, £12 = £21 = ^22 = 0, $i=0, $ 2 = <P 

K 1 it;{<f>})=[0 2 ~^ + Uoit)], K 2 (*;{*}) =2^. (53) 

Note, that the equation (^2|) determines a diffusion process for which <& (t) is continuous. 
Let the conditional probability meet the boundary condition 

P(*,t|*',t / ) = 5 ($-*')> (54) 



then for small time intervals it is easy to determine the solution of equation (52): 

1 



P (#,t|*', if) 



2vrVeAt 



cxp [*-*'- K ( <; <*» A< 1 T £ ^ (*-*'- K (*; {*}) Ai )} - 

2^S GXP ^ * - ^ -V 2 + Uo W) At ] 2 } ' * = *' + At ' 



(55) 



Thus, one can state that the evolution of the system in the functional space is characterized by a 

regular shift with the velocity K(t; {<&}), against the background of which the Gaussian fluctuations with 
diffusion matrix take place. As to the trajectory 3? (t) in the space P{*}, it is determined by the formula 
(see§) 
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* (t + At) = * (t) + K (t; {*}) At + F (t) At 1/2 . (56) 

The trajectory is seen from ( |56| ) to be continuous everywhere, i.e. <fr (t + At) $ (t), but is undifferen- 

tiable everywhere owing to the presence of term ~ At 1 / 2 . If we write the time interval in the form At = t/N, 
where N — > oo, then one can interpret the expression ( p5|) as the probability of transition from 3?^ = <fr (t') to 
= $ (t) during the time At in the model of Brownian motion. 



6 Solution of Fokker- Plank equation for the distribution of 6 coor- 
dinate in the limit t —> +00 



First, let us consider the equation fl2q). Taking into account ( |24| ) as well as the initial condition (23) one can 
have the solution of equation (EG) in the explicit form: 



<p = Q in exp [-2 J 9 (t') dt' I (57) 

— OO / 

Using (|57|) one can reduce the system of SDE (|25|)-(|2Q) to one nonlinear and nonlocal SDE of the form: 



' + 2 -O? exp 



1 (t') dt' +U (t)+F (t) = 0. 



(58) 



Using SDE fl5q ) one can obtain the equation for conditional probability Q(9,t) = P(0, t\0, 0) in the sense 
of distribution of coordinate 9 at the time t: 



d t Q(9,t)+d e J 



•{t')dt';t 



(59) 



with the flow of probability 



J 9, 



>(t')dt';t = - 



- fi£, exp -4 



'(t')dt' +C/b(t) 



Q - ed e Q 



(60) 



In view of the fact that in the equation ( |59| ) the factor of the shift [9 2 — £l 2 n exp (— 4x x J 9 (t') dt' ] + Uq (t) 

— oo J 

tends to infinity at the limits 9 — > ToOj a symmetric non-vanishing flows of probability J(— oo,+oo;t) = 
J (+oo, — oo; t) 7^ there arise on the borders. For solution of equation (|59|)-(|60|) one has to require the natural 
initial and boundary conditions 



lim Q (9; t) = S (9) , 



lim Q(6;t) = 0. 
\e\^oo 



(61) 



(62) 



The investigation of equation (|59])-(p0[) with initial and boundary conditions for arbitrary t is a very difficult 
problem. However, below we shall see that for the construction of averaged matrix in the model of roaming 
QRHO it is important to know the distribution of 9 coordinate when t — > +oo. 

Turning again to the solution of (|57]) note, that in the limit of large time intervals the self-averaging of 
trajectory 9 (t) takes place by virtue of the ergotic hypothesis, so that 



+oo 



lim if (t) = aQi 

t — > + OG 



a = exp 



' (t') dt' 



(63) 
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We should like to remind, that by virtue of the same ergodic hypothesis the constant a can be calculated 
by averaging over the ensemble. 

Now taking into account that for t — > +00 the density of flow reaches its limiting value 




J / = - lim J\0, 6(t')dt';t\ , (64) 



we obtain the equation for the probability of distribution of a stationary process, 

J of = [e 2 - a 2 fl 2 „ + U+] Q s - edeQ 8 , d e = d/ d6, (65) 



U+ = lim U (t) . 

t — >+oo 

In Eq.(|65|) the probability distribution for a stationary process is denoted by Q s that is obtained from 
Eq.@ to be § 

Q s (e,X, T 0) = e-^ 3 Q s (A, 7; 6) = 



e 2 ^ 3 Jo/exp 



dz exp [ ( — + A72: 



(66) 



where A = (fl m / e 1 ^) 2 , 7 = (flout/ '^m) 2 — a 2 + F / flf„ and = 6/ e 1 / 3 . The constant Jo/ is calculated 
from the condition of normalization of stationary distribution Q s to unity and has the form 



00 

J-/ = vr 1 / 2 ^ 1 / 3 J dz z- 1 ' 2 exp - A 7 z 



(67) 



One can obtain representation of Jo/ in terms of special functions. Passing to Fourier components in equation 
|6S|) we find M: 



j ; = [A1 2 (-a 7 ) + B1 2 (-a 7 )] , 

where Ai [x) and Bi (x) are linear independent solutions of Airy equation ]lT| 

// 

y - xy = 0. 



(68) 



(69) 



Note, that the normalization constant Jo/ has one important feature, namely, it defines the number of states 
with energy values E from —00 up to fl 2 ut per unit of distance. In the case when (i?/e 2 / 3 ) 3> 1 and E 1 > the 
expression (|67|) is calculated asymptotically: 



iV s = J 0/ P3 TT^E 1 ^ 



32 £ 3 



£ 6 



(70) 



One can use expression (|67| ) for determining of states distribution for concrete energy value. Taking into 
account the fact that bound states are characterized by negative energy values and making in equation ( |67|) 
the formal substitution fl 2 nf + Fq — * —E one can obtain expression for the number of states with concrete 
localization energy: 



\ L : n : -'/,/:: 1 ' 2 exp ( + 



(71) 



In the case when (E/e 2 / 3 ) ^> 1 and E > the expression ( |7l|) is calculated asymptotically and have the 
following form: 
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7r _1 -E 1/2 exp 



4E 3 / 2 
3 e 



1 + 



It is easy to calculate now the states distribution for concrete value of bound state energy: 



Pe 



Ne 



exp 



3 e 



1 + 



( £ 



\E 3 / 2 



(72) 



(73) 



It is worthwhile at the end of this section to pay attention to one important property of the stationary 
distribution Q Sl that follows, in particular, from numerical analyses. The fact is that for all 7 6 (—00, +00) the 
probability of positive values of the coordinate is higher than that for negative values (see FIG.l). It follows, 
hence, that a = and ip s = , where tp s is the imaginary part of complex coordinate $ in the stationary 
process limit. 



7 Calculation of the wave function of roaming QRHO 



Having at hand the expression for conditional probability (|55|) one can now pass to the averaging of the wave 
functional in the space, or, in other words, to the calculation of wave function of roaming QRHO. Writing 
the wave functional (|34|) in the moving reference system (see (|20|)), and taking into account the formula ^ for 
the average value of wave functional, we obtain 



*^(n|x, *;{*}) 



{*} 



- / D»{*}*<£>(n\x, <;{#}), 



(74) 



In this formula Dn{<fr} determines the total Fokker-Plank measure of the functional space Rs&x : 



Dfi{$} = dfi{& } ■ dfj,{& t } ■ lim 

iV— >OC 



J_ IN 

2ir V et 



N 



N 



(75) 



x n exp e-e'- (e 2 + ^ + u (t)) & 

fc=o L L 

where a, rf^t{#o} and d/i{<& t } are determined with the help of formulae 

a(e,t)=fDii{<*>}, 

dfJ. {*o} = 8 (Oo) S (Q in - ip ) d9 dip , 

dfi{^ t } = P(#,t| 0,0) rf^. 



d9k+idip k +i 



(76) 



In equations (|75|)-([76t) a(e, t) denotes the normalization factor for functional integral (74) with total Fokker- 
Plank measure and in the limit t — > +00 it acquires a constant value that, in general, is different from the unity. 
The integration over the measure dfi {$0} means the averaging in the initial distribution of complex coordinate 
$ at the moment of time t = —00. As for the integration over the measure dfi{<&t} , it provides the averaging 
in the distribution of complex coordinate $ at an arbitrary instant t . In particular, one can show that in the 
limit t — > +00 it assumes the following form: 



eW*oo} = d[i{& s } = Q s (e,X,'y;9)5((p)d9difi. 



(77) 



Now we shall try to study the (out) asymptotical state of the "oscillator+thermostat" system . 

Taking into account the relation (|^) one can separate the fluctuation process at the frequency (^) in the span 
of time (t c , +00), where t c is the zero time reference point in the (out) channel, by giving it with the help of 
white noise correlator with the diffusion constant e+. It is evident, hence, that the (out) asymptotic state of the 
"oscillator+thermostat" system will be also characterized by a complex probabilistic process ^ Q ut (m\x, t; {??+}) 
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in the extended space ^ out = Rl ut <8> R{ v+ }, where i] + = exp ifl out t+ J $ + (t')dt' is the corresponding 

RK ' V «« / 

solution of SDE (|9|). Proceeding with the same reasoning one can construct the wave function of the (out) 
asymptotic state 9br (m\x,t). 

At the end of this section we should like to note that in the limit e — + the measure of the functional 
space R{£} turns out to be of the Wiener Type, and the wave function of roaming QRHO, 9^ (n\x, t) steadily 
transforms to the exact solution <Ia+) (n\x,t) that describes the motion of regularly moving QRHO. 

8 The local stochastic matrix of transitions of roaming QRHO 

As was mentioned above, the roaming QRHO has two asymptotic states: 9 in (n\x, t) £ L 2 (Rj n ) when t — > — oo 
(see (0)) and respectively 9 out (m\x, t; {r) + }) G L 2 (Rl ut ® R(r} + )) when t — > +00. The goal of the theory of 
scattering is to construct a unitary operator that will change one of these states to the other pj . 

Prior to that let us consider the wave functional {m\x,t; {77}) that is the solution of SDE (|])-(0) and 
passes to the asymptotic wave functional 9 out (m\x, £+; {??+}) when t = t+ > t c , which in its turn passes to the 
wave function 9 out (m\x,t) when t = t c . Since the sets of functional 9^ (n\x,t;{£}) and 9^ (m\x, t; {rj}) 
form completely orthonormalized bases in L 2 , one can write the following decompositions: 

(n\x, t; {£}) =J2 S *» (*i ^ I*'! W}) (k\x, f ; W}) , (78) 

k 

(m\x, t; {„}) =J2 S *km (t'l U'} It; M) {k\x, t'; {£'}) , (79) 

k 

where the coefficients Sk n and St m are some random local evolutionary operators. Taking into account the 
orthogonality relations (E3) , one can obtain from (|7^ ) 



5 nm (i;U}|i^W}) = ^ ) W^*;m)*L ) (H^*';M)) ■ (so) 

One can easily find from expressions ( ^8[) and (|79|), with due regard for (^6|) and ([sb]), 

£ ^« (*; {CI I*'; {»/}) (*; iO \t'; W}) = <W, (si) 



k 

jstc 



that is basically a generalization of the unitarity condition for the stochastic matrix S s , the matrix elements 
of which are determined by means of the formula J8C|). 

Note that the unitarity of the stochastic matrix S stc (S stc ) + = I is a consequence of conservation of laws in 
the closed system " oscillator+thermostat" , that is defined in the extended space S. 

It is convenient to carry out explicit calculations of the elements of stochastic matrix ( |80| ) by means of the 
method of generating functionals. Let us consider the sum 



00 n 

(z\x, t; {£}) =E n=^c (n\x, i; {£}) , (82) 

n— 



where z is an auxiliary function. If we substitute the expression for wave functional ( |34j ) in ( jS2| ) and make 
an appropriate summation (see Jll|] ), then we shall have 

1 

*&> (z\x, t; {£}) = (^) 4 exp |-i (ax 2 - 2bx + c) J , (83) 
where the coefficients a, 6 and c are 
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&(*;{£}) = ^{wf) ' 

c(i;U}) = z 2 exp [-i2r{t-{W})\ 



(84) 



As it is seen from (|8^)-(|84|), the generating functional in x coordinate is a stochastic Gaussian packet that 
in the limit t — » — oo goes into an ordinary Gaussian packet: 



x exp — - ( Q,i n x z — 2^Qi n zxexp (—iQi n t) + ^ 2 exp {—i2Vti n t) + if^t 



(85) 



As for the generating functional of the (out) state, it is easily obtained from (|83|)-(84) by making simple 
formal substitutions z — > z + , n$„ — * f2 out , £ — > ?y + and t — > i + : 



* out (z + |x,< + ;{77 + }) = ( ^ ) ?y + V2 exp^ -- (a+a; 2 - 2&+x + c+) 



where the following notations are made: 



6 + (t +I z + ;{^ + }) = ^f& 



c+ (*+; {*7+}) = 4 ex P H2r+ (*+; i w +})\ 



(86) 



(87) 



r + ( t+iW )^ mt / M ^, }) . 
Now consider the following integral: 

J (z,t;{# 14, *+;{»?+}) = (*i+ ) (2|^t;{0)* o «t(^+N,i+;{r7+})) x . (88) 

Making appropriate substitutions for the generating functionals from (p3[)-(p^) and (p5[)-(p6[) and integrating 
over the coordinate x we obtain 



/ (z, t; {O K, i +; {,?;}) - (n in n „ t ) 1/4 

with the following notations: 

A(t;{£}|*+;{4}) 



U£4 



1/2 



exp 



.4 



(89) 



= -if +»%. 



S (t, z; U} |t + , z + ; {^}) + 



(90) 



C (t; {£} |t+; {v*+ }) = exp (-»2r) z 2 + exp (i2r + ) (z* + ) 2 . 
It is easy to show that I (z,i; {£} |z^_,t + ; is a generating functional for stochastic matrix elements: 



l(*,t;{£}l4>*+;{4})= E -^=5 nm (t; {?} {t?;}) 

— n V^™! 

7 1.771 — (J 



(91) 
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Making expansion in the left hand side of expression ( J9l| ) in the Taylor series in z and z* + we find 

1 



S nm (*;{£} l*+;{4}) = 



(*,*;{£} 14. *+; {4}) 



\Jn\m\ 

To calculate several first elements of the local stochastic transition matrix one can substitute (|9|) to (^2| 

1/2 

5 00 (t;{e}|t+;{4}) = ( 4 ^«n „ t ) 1/4 ex P (f){e4(er 1 -4(4)" 1 )} > 



(92) 



MM0K+;{4}) = [^oo (t;-[c> 1*+; {»?+})]» 

1/2 



S20 (*;{£}!*+; {4}) = 5 < 



on 



■ exp(-i2r)+(^) 



(93) 



S02 (t;U}|t + ;{4» =5, 



on 



ox|,(/2, + ) + (^) 1/2 e(4) ^oo 



One can see from expressions (|93|) that in analogy to the case of regular QRHO the transitions will occur 
only between asymptotic states of similar parity irrespective of the value of diffusion constant e. But in contrast 
to the regular case, it is characteristic that the symmetry in local matrix elements ( |93|) with respect to the 
transposition of quantum numbers of initial n and final m channels is violated. 

9 The averaged transition matrix of roaming QRHO. The probabil- 
ity of "vacuum- vacuum" transition 

The ultimate objective of our study is the calculation of transition probabilities of roaming QRHO. This aim 
in view it is necessary first to make the averaging of local stochastic elements (193]) . 
Definition 1. The expression 



nor 



lim 

t — ► +OO 
— > +OO 



S nm (t;{* + f}|t + ;{*++f + }) 



(94) 



{*}/ 



{*+} 



is termed the averaged (generalized) transition matrix of roaming QRHO, 

where the local S nm it; {3? + f} \t + ; {$ + + f + }) matrix element in ( p4| ) is obtained from ( p3| ) by writing the 
latter in an arbitrary movable reference frame with the help of transformation ( |20| ) . Remember that the first 
bracket in (p4) means a functional integration over the measure 



N 



Df (t) = lim II S (f (t fc ) - fo (t fc )) di (i fc ) , 



JV 



tk — M 



J k=0 



N ' 



(95) 



f (t) = {Re 



£o(t) 



; Im 



fo(t) 



The second bracket respectively means the integration over the measure 

N 

Di + (t) = Df+ (t) = lim J] 5 (/+ (i fc ) - / 0+ (ifc)) 4f+ . 

A*— °°fc=0 

fo+ = /o+ = {0; f^owt} ■ 



(96) 



In (Q) the third and forth brackets respectively stand for integration in and -R{* + } spaces. 

Note that the integration over variables f (t) and /+ (t) gives the average value of transition matrix in the 
steady (laboratory) frame of reference. It is also noteworthy that, as it follows from expressions (|8^)-(^l]), the 
generalized transition matrix S br is usually not a unitary matrix, i.e., S br (<S ,br ) + I. Nevertheless, there are 
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some limiting cases when one can simplify both the measures of functional spaces and R^ + j (by making 

these of the Wiener type), and the corresponding wave functionals. As a result of these simplifications, the 
contribution to the expression (Q) is made only by the integrals over the final distributions of coordinates 9 
and #_|_ and that, in its turn, makes unitary the generalized transition matrix Syy = I, where the index 

w means simplification of the matrix in the above sense. 

As an illustration of the proposed approach, consider now the probability of " vacuum- vacuum" transition 
in the case when the spaces and R(& + } have the Wiener-type measures. After the substitution of the 

expression for 5 o (t; {*} \t+; {*+}) to (|9j) and simple integration we find 



S b r (A, A+; p) = (1 - p) 1/4 {h (A, X+;p) - il 2 (A, A + ; p)} . 



where the following notations were made: 



(97) 



+ 00 



Ji,a(A,A + ;p)= / d9Q s (\, r ,6) / d9+Q s (X+, T ,h) 



a±l 
2a 2 



1/2 



(98) 



a(X;9\X+-9+) = 




1/2 



(99) 



In above formulae p stands for the reflection coefficient from a barrier in the one-dimensional problem of 
quantum mechanics with momentum K (x) — flo (x) , where t was replaced by x (see ) . As for the function 
j(p), for the frequency model in the form of step barrier 



O (i) = n in + {Q out - Clin) sign(t) 



(100) 



it has the form 



7(p) 



Q U 



.1/2 



1 



,1/2 



(101) 



Now one can write the final expression for the amplitude of " vacuum- vacuum" transition probability 



A^ (A, A+; p) = ^l^p { II ( A, A+ ; p) + 1% ( A, A+ ; p) } , 



(102) 



with due regard for the modification of the nature of fluctuation process in the "oscillator-thermostat" 
system. 

In case when e + — > 0, i.e., the final state of the system is described by the wave function, the expression 
(101 ) is strongly simplified 



AoLo (A; p) = V^~P {if (A; p) + If (A; p)} 



(103) 



h,2(X;p)=j dS (^r) Qs(\t,o), (104) 

— oo 

a(\,r,8) = U + ^ J ■ (105) 

The numerical calculations of " vacuum- vacuum" transition probability by the formulae ( |l03| )-( |l05| ) show its 
nonmonotonic behaviour (see FIG. 2). 
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10 Calculation of average values of dynamical variables. Thermo- 
dynamics of the nonrelativistic vacuum in asymptotic subspace 

r>l 

^as 

It is well known that the main object of quantum statistical mechanics is the density matrix that can be written 
in the nonstationary representation as (see fll9fl ) 

p(x,x',t) W m frn(x,t)ip m (x' ,t), (106) 

m 

where w m is the probability that at the moment of time t = the system is in the state ip m {x,t). The 
function (p n (x) is the solution of Schrodinger equation that satisfies the initial condition 

tp m (x) = ip m (x,t)\ t=0 , (107) 
At t = the following density matrix is defined with the help of a set of wave functions <p m (x) : 

p(x,x') =y^ w m (p m (x)<p m (x'). (108) 



Note that the density matrix in the form of ( |108[ ) was first defined by Dirac and von Neuman (see |2(J ) , the 
following microcanonical distribution being introduced for the coefficient w m : 

w m = exp f-^^j , (109) 

where E m is the energy of quantum level m, k is the Boltzman constant and T is the temperature of the 
thermostat. 

In particular, it follows from the definition ( |l06| )-( |l08|) that if at an initial moment of time t = the system 
is in the state ip m (x,0) with probability w m , then the probability for the system to be in the state tp m (x,t) 
at the moment of time t will be the same. One should note that the representation (108) is valid only if the 
quantum system weakly interacts with the thermostat. 

Below we shall study some relaxation processes that occur with the energy spectrum of vacuum-immersed 
quantum harmonic oscillator (QHO), as well as calculate the thermodynamic potentials of vacuum state when 
no any limitations are imposed on the amplitude of interaction between the quantum oscillator and vacuum. 

Because each quantum state m in the problem under consideration is described by a complex random process, 
it makes sense to apply here the thermodynamic description. 

Definition 2. The partial density matrix is defined to be 

piTj = (x, t; {£} \x', t>- {£'}) - {* stc (m\x, i; {£}) * stc (m\x',t' ;{?})} . (110) 
Let us remember, that the wave functional ^ s tc (m\x, t; {£}) in ( |110| ) describes QHO state with the frequency 

O(t) = n as + n+{t;{W+}), Q as = const, (111) 
where the stochastic term of frequency squared have the following form 

U+ (t) = 2Cl as Q + (i; {W+}) + Q 2 + (t; {W+}) = F 0+ + F+ (t) , (U+ (t)) = F 0+ = const . (112) 
We shall suppose that 

(F+(t)) = 0, (F + (t)F + (t'))=2e + S(t-t / ). (113) 

Definition 3. The mathematical expectation of the stochastic operator A (x,t; {9}) in the quantum state 
m is defined as 

A= t^+o { Spx Sp {e}^*te I s Vx Sp {i} p^ t n J j. (ii4) 
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where Sps x \ means the integration along the diagonal line x = x', and Sp^gy denotes respectively the 
functional integration along the diagonal {9} = {9'} with measure (f75|)-([77|) taking into account that in the 
asymptotic subspace R as the following notations were made fio(i) - * ^as, Fq — > i*b+, 7 7+ = 1 + Fo+/ ^Li 
and e — > e+ . 

Definition 4. We shall call the following distribution function "the nonequilibrium partial" one 



0(e+,E,t) = S Px {Spwp™} = S Px {s P{e} p ( ™J (x,t;{e} \x', f ; {9'})} , 



(115) 



where E is the average value of energy in the quantum state characterized by the index m, Sp^gj standing 
for the functional integration along the diagonal line {9} = {9'} but with the measure 



Dfi {9} = a- 1 ( e+ , n as ; t) ^f±lM Dfx {6} , 



(116) 



where Ne(s+, 9; t) is the solution of the Fokker-Plank equation (p9|)-(|60|) after replacement of Uo(t) by —E, 
E > 0, and iVs, respectively, after replacement of Uo(t) by E. It is easy to see that the ratio Ne/ N% gives the 
distribution of quantum states in the vicinity of binding energy in the function al sp ace R{g + y. It is necessary to 



stress that the average energy E here is determi ned w ith the help of formula 



114]) and, hence, depends on the 



quantum number m. After simple integration in ( 115 ) with due regard for (116) and passing to a limit t 
we obtain for the equilibrium distribution function of quantum states 



-co 



^Hs + ,E) = ^ 



E>Q, 



Ny, = N- 



(117) 



where 



dzz- 1 ' 2 



exp 



z 
12 



E 

JV2" 



(118) 



Having in view the asymptotic estimates (|70|)-(|73|) one can conclude that when the condition Ej e 2 / 3 <C 1 is 
met, the partial distribution function z?( m )(e + , E) passes into the microcanonical distribution when the constant 
e + = jE^kT. Owi ng t o that, the integration of stochastic density matrix ( 110 ) over the functional space R{g + } 
with the measure (116) and summation over the index m gives the density matrix in the representation of Dirac 
and von Neuman (106)-(10£). 

Now, having the partial distribution function i?^ m '(e+, E) one can determine the thermodynamic potentials 
of the specific level m : 

a) the average internal energy 



U (ra) (e+,E) = -d e+ {ln^ (£+,£)} 



b) the free Helmholtzian energy 



F< m ) (e+, E) = -E+ 1 LatfM (e+, E) , 



c) the entropy 



S( m » (e + , E) = £ + fc{uW (e+,E) + F'™' (£+,£)} 



(119) 



(120) 



(121) 



Having in view the illustration of results we shall carry out the calculations for vacuum state, i.e., m = 0. 
An appropriate partial stochastic matrix in the immobile reference will have the form (see Section 4) 



p ( 2 (x, t- {9} | x', f ; {9'}) = (IW) S exp |_^ (x 2 + x ny 
\ } 9{t)cLt-\ J 9{r)dT - i [9(t)x 2 - 9{t')x' 2 } 1 , 



(122) 



where fl as is the permanent frequency of QHO in the asymptotic space Rg 8 . 
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Substituting the expression for stochastic potential (|^) and stochastic density matrix ( |122| ) to flll^ ) and 
making simple calculations we obtain the following expression for the average value of energy of "vacuum- 
oscillator" system in the ground state: 



(A+;fi„.) = \Sl as J+ J d^- 3 / 2 exp f-^-X+t 
o v 

+ ±f! as |l - i J dzz3/2 exp (-g - A+z) | + 

+i V>7 J + ^ dzzV2exp (-12- A + Z J ' 



(123) 



where the notation was made 



JI 1 = tt 1 /^! 173 / dzz-^exp 



A 4. z 

12 + 



Note that the expression (123) was obtained taking into account that (F)^ w y = 0. One can see that the 
first term in (123) is diverging that corresponds to the infinite energy of physical vacuum. The second term 
corresponds to the energy of oscillator in the ground state that is shifted as a result of interaction with the 
nonrelativistic vacuum (see FIG. 3) 



Sl(A + ;fi os ) = ±n as {l-£ \(d a \nA(-X + + a)) 2 + d 2 a lnA(-X + +a)] } 
A (-A+ + a) = Ai 2 (-A+ + a) + Bi 2 (-A+ + a) , 



(124) 



It is noteworthy that the second term in (124) is an analogue of the Lamb shift of the energy level, that is 
well known from the QED |2l|]. The third term in ( jl23| ) corresponds to a broadening of the ground level and, 
hence, is inverse proportional to its decay time 



{d^ln/lt-A + a)}- 1 



a=0 



(125) 



Now, based on the formulae ( 1 17 )-( 121 ) , one can write an expression for the entropy of QHO in the m-th 
quantum state immersed in the physical vacuum 



S^l {e + ,E 0SC ) = -ke + \d a \n 



A 



2/3 

-+ 



A 



JSJW 



-e- 1 In 



A 



_2/3 



A 



E osc 
2/3 



(126) 



To obtain an expression for the entropy of oscillator in the ground state one has to make a substitution 
E osc — > E^osc in the formula (126) 



(X 0) ) 



k{0. 



(0) 



il/2 



(-/*£>) 



il/2 



^-i/2(-/3. 



(on 



hi 



^- 1 /2(^° ) ) 

^4-1/2 (-M ' 



where the notations were made 



00 



(o) _ 



2/3 • 



(127) 



(128) 
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11 Conclusions 



At present the three following major schemes for the quantum chaos initiation are discussed: 

a) the dynamical one, when the classical analogue of the quantum object under study is a nonintegrable 
system ||. An evident example of such a case is the three-body problem |l^|-||T5|; 

b) the measurement problem comprising the issues of mesophysics, of the irreversibility problem and quantum 
jumps |, 

c) randomness, and, therefore, the irreversibility as the basic properties of the physical world, owing to which 
the really existing object,-the physical vacuum (see |^]), is included in theoretical consideration. 

The objective of the present paper was to investigate the third case. The main point here is the idea that 
the quantum object and the physical vacuum are considered to make a joint system described in the framework 
of single equation. Since such a system has infinite number of the degrees of freedom, it is convenient to 
mathematically formulate the problem in the framework of complex SDE in the extended space S = R 1 (g> R{^} ■ 
In this case the system "quantum object+thcrmostat" will be described by the wave functional, - a complex 
probabilistic process, and not by the wave function. As a specific example, the problem of one-dimensional 
QRHO motion in the Euclidean space R 1 has been considered in the work. It was shown that in the complex 
SDE (|lj)-(Q) the variables were separated by means of Langevin-type standard real nonlinear SDE and the 
solutions for the wave state were obtained in the form of complex orthonormalized functionals (n \ x, t; {£}) 
in the space L 2 (R 1 F° r nonlinear Langevin SDE the corresponding Fokker-Plank equation has been 

obtained and with its help a complete positive measure for the functional space was constructed. This last fact 
permits the construction of uniformely converging representation for the averaged wave function [n \ x, t) of 
moving QRHO on the conventional space-time continuum (R , t) by means of functional integration of the wave 
process \&g£? (n | x,t;{£}) in the space R{£}- The obtained mixed continual-undulatory representation unites 
two concepts that at the first sight seem contradictory: the quantum analogue of the Arnold transformations 
[ fil| that excludes the rise of chaos inside the given trajectory beam, i.e., forbids any changes of beam topology, 
and the functional integration that allows for the contribution of various topological pipes and thereby assists to 
the generation of the chaos. With the help of wave functionals ty^c {n \ x, t; {£}) an expression for the elements 
of transition matrices (see (|30|)) was constructed. It was shown that the stochastic matrix formed in such a 
way is unitary S stc (S stc ) + = I. The generalized transition matrix S br for QRHO was obtained by averaging 
of the stochastic matrix S stc in the R{£} space and it was shown to be generally nonunitary. In other words, 
in the theoretical scheme in question the quantum oscillator is an open system in the space-time continuum 
(i? 1 ,^) and here the conservation laws are violated. Nevertheless, when the oscillator weakly interacts with the 
thermostat the transition matrix is simplified and turns to be unitary, S^(S^) + — I, i.e., the conservation laws 
are effective again. 

In the present paper a detailed analysis of the probability of " vacuum- vacuum" transition AqI^q (A, p) is 
given in terms of A that characterizes random fluctuations of the thermostat before and after the crossover, 
as well as in terms of barrier reflection coefficient p in the relevant problem of the one-dimensional quantum 
mechanics. It was shown that for relatively strong interactions of the oscillator with the thermostat, i.e., for 
relatively small parameter A, the p-dependence of the probability is nonmonotonic (see FIG. 2). 

Based on the example of vacuum-immersed QHO, a novel f orma lism of quantum statistical mechanics has 
been developed in the framework of density matrix approach (|llO| ). It was shown that in the limit of weak 
oscillator- vacuum interaction, i.e., at A + — > oo it can be reduced to the Dirac and von Neuman density matrix. 
In the framework of novel approach the relaxation effects related to the energy spectrum of the oscillator were 
studied in detail. In particular, analytical expressions for the energy of ground state with broadening and shift 
(analogous to the well known Lamb shift from the Quantum Electrodynamics) as a function of constant A + , 
and for the entropy of the ground (vacuum) state of the oscillator were obtained. 

We should like to draw attention to one very important property of the quantum representation under 
consideration, viz., that it permits to relate the domain of quantum chaos with that of classical chaos ||l5| , the 
transition to the classical chaotic dynamics taking place in the limit h — * when e =/= 0. In case when e — > 
and h ^ SDE (Q) goes into the usual Schrodinger equation and describes a reversible quantum process. The 
classical theory of reversible processes is obtained when one tends to zero Ti — > and e -— > 0. 

An explicit relation between the quantum mechanics and quantum statistical thermodynamics was also 
established to specify, in particular, the range of applicability of the statistical matrix of Dirac and von Neuman. 

In conclusion we should like to note that one can successfully develop the quantum mechanics in the frame- 
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work of complex SDE (|IJ) both for the model of multidimensional oscillator and also for other exactly solved 
nonstationary problems of quantum mechanics. Its covariant generalization is also not difficult, and we postpone 
the detailed discussions of these problems to our later publications. 
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Figure 1: The dependence of the distribution of stationary process Q s (A; 0) of over 9 for different values of 
parameter 1/A <~ e. 
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Figure 2: " Vacuum- vacuum" transition probability in dependence of A and p. 
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Figure 3: Dependence of oscillator ground state energy, its shift and entropy of vacuum over the parameter A + . 
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